Cavity QED simulation of qubit-oscillator dynamics in the ultrastrong coupling regime 
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We propose a quantum simulation of a two-level atom coupled to a single mode of the electromag- 
netic field in the ultrastrong coupling regime based upon resonant Raman transitions in an atom 
interacting with a high finesse optical cavity mode. We show by numerical simulation the possibility 
of realizing the scheme with a single rubidium atom, in which two hyperfine ground states make up 
the effective two-level system, and for cavity QED parameters that should be achievable with, for 
example, microtoroidal whispering-gallery-mode resonators. Our model takes account of dissipation, 
and we pay particular attention to observables that would be readily observable in the output from 
the system. 

PACS numbers: 42.50.Pq,05.70.Fh,03.67.Ac,42.50.Ct,32.80.-t,42.60.Da 



I. INTRODUCTION 

The idea of analog quantum simulations, that is, engi- 
neering one quantum system to replicate the behaviour 
of another that one wishes to study, was first proposed 
by Feynman [l| with regards to the infeasibility of sim- 
ulating quantum systems on classical computers. The 
simulating system typically involves a larger number of 
degrees of freedom, and the effective simulation relies on 
precise variation of system parameters through exquisite 
experimental control. This idea has become reality in 
recent years 0, Q and fundamental models of interact- 
ing quantum systems have been realized thanks, for ex- 
ample, to advances in the control and manipulation of 
ultracold gases in optical lattices Systems based on 
cavity quantum electrodynamics (cavity QED) also offer 
exciting possibilities, as exemplified by the recent demon- 
stration of the Dicke quantum phase transition with a 
superfluid atomic gas in an optical cavity [11 , in which a 
pair of discrete momentum states of the gas was used to 
simulate the two-level atoms of the original Dicke model 
of superradiance 0, 0] ■ 

In this work we focus on a closely related model - the 
Rabi model ~ which describes the interaction of a two- 
state system (atom or qubit) with a single quantized har- 
monic oscillator (electromagnetic field mode) through the 
Hamiltonian 

Hi^ = -^a^_ + Loa) a + g{(j+ + (j^){a + a^), (1) 

where wq and ui are the qubit and oscillator frequencies, 
respectively, g is the interaction strength, {(Jz,(J±\ are 
two-state operators, and a (a^) is the annihilation (cre- 



ation) operator for the oscillator. This Hamiltonian pre- 
dicts accurately many physical situations where an atom 
- artificial or real - is interacting with a confined cavity 
field. However, when this model is used in the descrip- 
tion of conventional cavity QED systems, it is generally 
the case that the coupling constant g is much (i.e., or- 
ders of magnitude) smaller than the frequencies loq and 
Lo. This means that terms in that do not conserve 
the total excitation number can be neglected, in what is 
known as the "rotating wave approximation." This leads 
to probably the most studied model in quantum optics, 
the Jaynes-Cummings model Q: 
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Recently, however, going beyond this approximation has 
gained newfound interest, as novel experimental systems 
have pushed the boundaries for the strength of the cou- 
pling between (artificial) atoms and cavity field modes. 
In particular, in experiments using circuit QED Q and 
semiconductor microcavities, couplings between artificial 
atoms and cavity modes have reached the so-called "ul- 
trastrong" regime [l0l - [l3j . Here, values of g/uj ^ 10% 
have been realized, so that g is no longer small compared 
to the mode frequency. 

These experimental advances have in turn stimulated 
new theoretical investigations (see, for example, (T^ - [l9j ). 
as the conventional treatment is no longer adequate. 
To appreciate this necessity, notice that for the Jaynes- 
Cummings model, Hjc conserves the total excitation 
number (i.e., \{(Jz + 1) + a) a), and, together with conser- 
vation of energy, this means that the spectrum is easily 
found. If, on the other hand, the terms aa- and a^cr+ 
cannot be ignored, then this is no longer the case, and 
until recently no analytical solution was known. This was 
remedied in [l^ by Braak, who showed that a discrete 
symmetry, the conservation of parity H = — cr^(— 1)*^ '^j is 
enough for the model to be integrable and derived ana- 
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lytical expressions for the spectrum (see also [l9|). Even 
though some intuition can be drawn from the conserva- 
tion of parity [13, HH , the conventional picture of excita- 
tion exchange between the atom and the field is clearly 
not valid, and the spectrum of the Rabi model and the 
associated dynamics has been shown to exhibit a vari- 
ety of novel and significant nonclassical effects, such as 
initial state revivals [20,, j strong atom- field entangle- 
ment (23| . and generation of photons from the vacuum 
through modulation of the coupling constant Q or qubit 
frequency [25l[26|. 

Based on the above discussion, it follows that practical 
realizations of the Rabi model in the various regimes of 
coupling strength will be important to completing our un- 
derstanding of the coupling of light to matter at the most 
fundamental level. We suggest that quantum simulation 
can be a valuable tool in this respect. Indeed, an analog 
simulation of the Rabi model based on light transport in 
waveguides was realized in |27| , and a way of simulating 
the model based on circuit QED was recently proposed 
in psj . Such effective realizations are important not only 
because they have the potential to simulate even greater 
coupling strengths than can be achieved at present by 
"direct" coupling between two-level systems and cavity 
modes, but also because they offer flexible means of con- 
trolling other model parameters, as well as convenient 
ways of probing the dynamics through well-defined out- 
put channels and measurements. 

In this vein, we wish to propose a realization of the 
Rabi model with a single real atom coupled to a high fi- 
nesse optical cavity mode. Two stable hyperfine ground 
states of the multilevel atom make up the effective two- 
level system, while (resonant) Raman transitions between 
them are induced by the cavity field and auxiliary laser 
fields. This allows us to simulate the Rabi model with 
essentially arbitrary tuning of the effective frequencies 
and coupling constant, so that any regime, including ul- 
trastrong and deep-strong coupling, can be accessed and 
explored. Although the model could in principle be real- 
ized with a variety of (alkali) atoms, we focus on the Dl 
line of ®^Rb, and consider cavity QED parameters that 
should be achievable with, for example, microtoroidal 
whispering-gallery-mode (WGM) resonators [l^. 

Inevitably, all such systems are subject to dissipation, 
and the Hamiltonian model must be expanded to an open 
systems treatment. We will see that our effective model 
is described by a master equation of the form 

p = —i[Hfi, p] + K{2apa) — a)ap — pa^a), (3) 

where k is the cavity decay rate. On the one hand, we can 
realistically expect k to be small enough for the Hamil- 
tonian dynamics to dominate on appreciable time scales, 
so that the Rabi model dynamics is prominent. On the 
other hand, the dissipative cavity QED setup provides 
us with a convenient means of observing the dynamics 
via the output field of the cavity. Furthermore, the open 
systems dynamics is of interest in itself, and the steady 
state of Eq. ([3]) is known to possess some very interesting 



features [30|. It should be noted, however, that such a 
master equation is incorrect for circuit QED s yste ms gen- 
uinely in the ultrastrong coupling regime mM 111,131, 

as it predicts unphysical results such as photon genera- 
tion from the vacuum. In our case however, such a pho- 
ton fiux is perfectly legitimate, given the energy input 
from the laser fields that help to drive the Raman tran- 
sitions that we use to implement our simulation of the 
Rabi model. Our scheme thus offers the possibility of 
realizing an open systems version of the Rabi model that 
is subject to a simple ("conventional") interaction with 
the environment (i.e., cavity damping) and can hence 
explore certain novel features of the Rabi model beyond 
what might be achievable with genuine systems. 

In addition, our scheme offers a generalization of the 
Rabi model by the addition of an extra term in the ef- 
fective Hamiltonian describing a non-linear coupling be- 
tween the qubit and the oscillator. This term takes the 
form 



(4) 



where gz is a tunable parameter, and can be viewed 
as a dynamical shift of the oscillator frequency, uj — >■ 
Lj + gzCz/^. A coupling of this nature typically arises 
as an approximation to the Jaynes-Cummings model in 
the dispersive limit (i.e., when jcij — wqI ^ g), whereas for 
our model it can be present independent of the relative 
sizes oi ujQ, w and g, and of a magnitude that may in 
fact be comparable to or even larger than these param- 
eters. Although our initial focus will be on the effective 
realization of the Rabi model, Eq. ([1]), such that we set 
gz to zero or to a small enough value the effect of the 
non-linear coupling is unimportant, we will see that vari- 
ation of this parameter is extremely interesting in itself 
and can lead to dramatic changes in the system's prop- 
erties. In particular, we find two "critical" values of gz 
{g'l ~ ±2a;), about which, in the deep strong coupling 
regime {g/ui > 1), sharp transitions occur in the quan- 
tum state of the system, marked by clear signatures in 
the cavity field. 

Finally, we would like to comment on the connection 
between the model proposed here, and models simulat- 
ing the collective interaction of an ensemble of two-level 
atoms with a single field mode. Our scheme is essentially 
the single atom version of the proposed realization of the 
Dicke model presented in lH . By showing that an effec- 
tive ultrastrong coupling regime is in principle achievable 
on the single atom level using a '^^Rb coupled to a high 
finesse optical microcavity, we also offer a specific means 
of realizing the many-atom Dicke model quantum phase 
transition as proposed in (33| (which requires only strong 
collective coupling of TV atoms to a cavity mode, a much 
weaker requirement than the very strong single-atom cou- 
pling assumed in this work). 

In fact, the above-mentioned experimental realization 
of the Dicke quantum phase transiti on pi was carried out 
using a scheme analogous to that of [33| , but based upon 
laser-plus-cavity-mediated, resonant Raman transitions 
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between discrete momentum states of a Bose-Einstein 
condensate. Notably, that particular scheme also gives 
rise to a nonlinear coupling term of the form J^a^a, 
where Jz is the many-atom inversion operator. The Dicke 
model including this term has been studied theoretically 
in the thermodynamic limit (TV — > oo) in (s^. [35l|. where 
phase diagrams for the scmiclassical steady state have 
been mapped out for parameters related primarily to the 
experiment of [Ej. The nonlinear atom-photon coupling 
was shown to be of fundamental importance; in particu- 
lar, a new superradiant phase is possible if the effective 
nonlinear coupling constant (g^) is negative. Our scheme 
could also offer a flexible platform for exploring such a 
regime and, indeed, the sharp transitions we already see 
at the single-atom level highlight this possibility. 

II. THE MODEL 

A. Full system 

The physical configuration that we consider here cm- 
ploys electric dipole transitions on the Dl line of a single 
^''Rb atom. By coupling the atom simultaneously to an 
optical cavity mode and two laser fields, two stable hy- 
perfine ground states - one in the F = 1 level and one 
in the F = 2 level of the 5^5*1/2 state - are connected 
through a pair of (distinct) Raman transitions. The spe- 



cific scheme is illustrated in Fig. [TJ In particular, err- 
and (T_ -polarized laser fields, separated in frequency by 
approximately twice the ground-state hyperfine splitting 
of 27r • 6.835 GHz and each far from resonance with the 
— 5^Pi/2 transition frequency, combine with a tt- 
polarized cavity mode to drive resonant or near-resonant 
Raman transitions between pairs of states, each of which 
consist of one state in the F = 1 level and one in the 
F ~ 2 level. We will eventually focus on just one pair of 
states, \F = 2,m = —2) and |F = l,m = —1), but for 
the moment consider the most general model. 

We introduce atomic dipole transition operators A''p\,, 
connecting level F in the 5^5*1/2 state to level F' in the 
5^Pi/2 state with polarization p: 

F 

A%= {F.m\l^p\F\m+p)\F,m){F\m+p\. 

m — — F 

(5) 

Here, m labels the magnetic sublevel, p = {— 1,0, +1} 
denotes {(T_, tt, cr+}-polarization, respectively, and 
/ip is the corresponding dipole operator. The 
dipole matrix elements are normalized such that 
^Pp\{F,m — p\iJLp\F' ,ni)\'^ ~ !> ^^'^ numerical values 
for the elements can be found, for example, in [3^. We 
can then compactly write the full Hamiltonian of the 
system in the form (setting h — 1) 



Hdi —Hq + H', 



Ho=uJr.a^a + uj2 ^ |2, m) (2, m| + ^ to) (l', m| + cj^ ^ |2', m) (2', to|, 



H' 



9JA'l^>+A-;,>+A^^>+A 



22 



f}+H. 



(6) 
(7) 

(8) 



r 



Here, uJc is the cavity frequency, uJi (i = 1,2) is the fre- 
quency (energy) of the atomic F = i level relative to the 
{^^Si/2,F = 1} level, with primed quantities denoting 
the excited levels, and ujli and uj^^ are the laser fre- 
quencies. For brevity, we use the notation = |_F = 
i,TO = n), \i' = \F' = i,TO = n), while a (a^) is 
the annihilation (creation) operator for the cavity mode. 
The laser Rabi frequencies are fii and f22, while gc is the 
atom-cavity coupling strength. Note that the terms in 
H' that arc proportional to {A^~^^\ A^j^^\ A^y^^\ A^^^'^} 
(and their conjugates) do not participate in any resonant 
or near-resonant Raman transitions and arc not shown 
in Fig. (TJ but these "off-resonant" terms can induce non- 
negligible shifts of the effective two-level frequency split- 
ting that we derive below. 



Including cavity (field) decay and atomic spontaneous 
emission at rates k and 7, respectively, the evolution of 
the system density operator, p, is given by a master equa- 
tion of the form 



p,F,F' 



A 



(p) 

FF' 



(9) 



where V[0]p = 20pO'f - O'^Op - pO^'O for any opera- 
tor O. The spontaneous emission rate for the Dl line of 
^^Rb is 7/(27r) = 5.7 MHz. We are most interested in 
parameter regimes where the atomic excited state pop- 
ulations are negligible, so we will neglect atomic spon- 
taneous emission in the effective model that we derive 
below. However, spontaneous emission will be included 
in all of our numerical simulations of the full model above 
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and we consider its effects in section Hill 



B. Reduced system 



As alluded to above, we assume very large detunings 
of the fields from the atomic excited states, so that these 
states are only ever virtually populated (assuming no ini- 
tial population) and can be adiabatically eliminated to 
yield an effective model involving only ground states (i.e., 
levels in the state 5^51/2)- To do this it is first conve- 
nient to move to a rotating frame through the unitary 
transformation defined by U{t) ~ exp(— iiJoi) with 



Ho =(i^L2 + 1^2)0^ a + (1)2 ^ |2, to) (2, to| 

711 

+ C^i,^|2',TO)(2',TO|, (10) 
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FIG. 1: (Color online), (a) Energy level diagram for the Dl 
transitions (not to scale). The two laser fields (blue dashed 
lines) are (7_- and o-+ -polarized, respectively. The cavity 
mode (red solid lines) is vr-polarized and participates in res- 
onant (or near-resonant) Raman transitions from both the 
F = 1 and F — 2 ground states, (b) Detailed view of the 
leftmost part of the energy level diagram. Horizontal dashed 
lines indicate detunings of the fields from the excited states. 



where (I>2 = {uili —^l^)/'^ is a frequency close (or equal) 
to cl>2, the ground state hyperfine splitting. We will from 
now on assume this transformation when referring to 
Eq. (O. Defining 



Ai =^2 - WLi, 

A2 - (WL2 +^^2), 



(11) 

(12) 



the condition for the validity of the adiabatic elimination 
is that IA1.2I ^ |^^iM^^2|, 't, 7- Neglecting sponta- 
neous emission and terms rotating (in the transformed 
frame) at frequency W2, we obtain an effective Hamil- 
tonian describing energy shifts of the ground state levels 
due to the various fields, and (Raman) couplings between 
pairs of levels, each pair consisting of one level in the 
F = 1 state and one in the F = 2 state. Focusing on the 
two leftmost levels as depicted in Fig. [IJb) , the relevant 
part of the effective Hamiltonian can be written 



^^off ~^ '^'^^^ + {Qi^+o- + .92C-a + H.c.) 

H CTzO' a. 

2 



(13) 



Here we have identified \g) = |2, —2), |e) = |1, —1), intro- 
duced Pauli operators a, = \e){e\ — \g){g\, (j+ = {cr^y = 
\e){g\, and dropped constant energy terms. The param- 
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eters of the effective Hamiltonian are given by 

Wo 



2 Ai 2 A2 6A2+w^i 

12 A2 + W2 12 Ai + 2w2 ^ ' 

/I 1 1 1 1 1 \ 

+ + (15) 



1 
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32 
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Ai ' 



3. =l3c|^ 7 



1 



1 1 



1 1 



where 



and 



4 A2 + w^i 12 A2 3 Ai 

Sc ^OJc - (wLa + ^^2), 



(16) 
(17) 
(18) 



(19) 
(20) 



— 0J2 — oj'i = 27r • 812 MHz is the excited state 
hyperfine sphtting. Note that, in fact, A2 + 1021 = Ai + 
u]2- The numerical prefactors to the various terms are 
products of dipole matrix elements associated with the 
atomic transitions 1361. 



C. Realization of the Rabi model 

We can choose gi = g2 = g (real) , so that H^s reduces 
to the form of a generalised Rabi model. 



HcS = — (Tz + OJO^a + g{a+ + CT_)(a + a)) 



9z t 



(21) 



The term proportional to gz can be made small com- 
pared to the other terms (or even zero) with a judicious 
choice of the physical parameters. Provided this is so, 
Hgg achieves an essentially faithful realization of the Rabi 
model. Furthermore, the effective frequencies and cou- 
pling strength of the model are all determined by either 
level shifts or Raman transition rates, which are tunable 
via the laser frequencies and intensities and can therefore 
be chosen to be of the same magnitude. In other words, 
we have a model of a two level system coupled to a single 
mode of the electromagnetic field where it is in princi- 
ple possible to access any regime of coupling strength - 
strong, ultrastrong or deep-strong coupling - as defined 
by the ratios ujq/uj and .g/w. 

Of course, our realization is with an open system and, 
including cavity dissipation, the master equation for the 
evolution of the reduced system density operator is 



(22) 



For the full Rabi model dynamics to be observable, we 
clearly require that the parameters of the model exceed 



the cavity field decay rate k. Hence, our proposed real- 
ization demands a strong-coupling cavity QED system, 
i.e., gc ^ K, so that, in particular, the effective coupling 
strength, g, can also exceed k. 

Cavity QED with microtoroidal optical resonators, in 
which atoms close to the surface of a microtoroid couple 
to the evanescent fields of whispering gallery modes (see, 
for example, [l^,!!^), should be capable of providing such 
a system. Very small mode volumes offer the prospect of 
electric dipole coupling strengths reaching values in the 
hundreds of MHz, while ultrahigh quality factors exceed- 
ing 10^ correspond to sub-MHz field mode decay rates at 
alkali atom transition frequencies (2^. Below, we choose 
the value gd (27r) = 200 MHz for most of our numerical 
simulations of the full model, and consider field decay 
rates in the range k/[2tt) = 0.01 - 0.2 MHz. 



D. Variation of Rabi model parameters 

To demonstrate how the effective parameters, 
Eqs. p4)) "([T8 )) . can be tuned through a variation of the 
physical parameters, we consider a situation where the 
cavity coupling constant gc is fixed. The effective non- 
linear coupling gz can then be set by an appropriate 
choice of detunings Ai and A2. The linear coupling 
strength, is tuned, for example, by varying fii, with 
^2 chosen so that gi = g2, which amounts to the con- 
dition ^2/^2 + ^2/{^2 + ^21) = 217i/Ai. This also 
sets the values of the quantities — <5 and uj — 5c- 
The effective frequencies and cj can then be tuned 
through variations of the parameters 5 and 5c- Since 
{\5\, \5c\} ^ {|Ai|, |A2|,a)2}, the adjustment of cjo and u 
in this manner is essentially independent of the tuning of 
g and gz- Note also that 5c could be adjusted simply via 
the cavity frequency Uc, while an external magnetic field 
could also control W2 and hence 5 through the relative 
Zeeman shift of the levels |2, —2) and |1, —1). 

Variation of the effective parameters as a function of 
laser Rabi frequencies and detunings is illustrated in 
Fig. [2] for gc/(27r) = 200 MHz. For simphcity, we set 
5 = 5c = Q, but note that the curves for cj and loq are sim- 
ply translated (uniformly) up or down for finite values of 
5 and 5c- In the top panel, the detunings Ai/ (27r) = —26 
GHz and A2/(27r) = —20 GHz are fixed (which sets a 
value for 5^/(27r) = -0.18 MHz), while the Rabi fre- 
quency r^i (and D.2) is varied. In the bottom panel 
Q,i/{2tt) = -320 MHz is fixed, while A2 is varied; VI2 
and Ai are varied correspondingly to satisfy gi = g2 and 
A2 + UJ21 = Ai + a;2, respectively. 



III. THE RABI MODEL SIMULATION 

In this section we first investigate the validity of the 
effective model, Eq. (j22p . by comparing its numerical so- 
lution to that of Eq. ^ with suitably chosen parame- 
ters [3^ . We focus on the simulation of the Rabi model. 
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FIG. 2: (Color online). Effective parameters in units of 
(27r)-MHz: g (blue, solid line), cjo (blue, dashed line), lo 
(black, dash-dotted line) and (black, dotted line). Top 
panel: Ai/(27r) = -26 GHz, A2/(27r) = -20 GHz, while 
fii is varied (VLi is varied simultaneously in such a way that 
g\ = gi for all fii). Bottom panel: r2i/(27r) = 320 MHz, 
while A2 is varied, with Q.2 chosen such that g\ = 172 for all 
A2. For both panels gc/i^n) = 200 MHz and 5 = 5c = (so 
A2 — Ai = a;2 — a;2i). 



Eq. ([T]) (i.e., we choose small values of gz - larger val- 
ues will be considered in Section |IV| . in the ultrastrong 
{g/uJ ^ 1) and deep strong couphng (g/cu > 1) regimes, 
with a focus on the behavior of the mean intracavity pho- 
ton number and the initial state revival probability. Hav- 
ing established a valid operating regime for the effective 
model, we then investigate briefly field-atom entangle- 
ment in the (dissipative) Rabi model. 

The open system nature of our setup makes it espe- 
cially accessible to experimental investigation. We can 
imagine, for example, a situation in which the atom is 
prepared in one of the states \g) or |e), while the cavity 
mode is initially in the vacuum state. Then, the lasers 
are turned on so that the system evolves according to 
Eq. ([9]). The cavity output can be continuously moni- 
tored by photon detectors, such that one can infer, for 
example, the intracavity photon number, field quadra- 
ture amplitudes, or photon correlations. The interaction 
can be stopped at any time by simply turning off the 
laser fields, after which one could also measure the state 
of the atom by, e.g., fluorescence detection. 



Cavity and laser parameters 



Set 


9c 




^2 


Ai 


A2 


S 


5c 


la 


200 


-160 


-120 


-26000 


-20000 


0.19 


0.40 


Ha 


200 


-320 


-240 


-26000 


-20000 


0.77 


0.40 


Ilia 200 


-640 


-480 


-26000 


-20000 


3.1 


0.40 


lb 


200 


-100 


-65.0 


-17000 


-11000 


0.25 


-0.034 


lib 


200 


-210 


-130 


-17000 


-11000 


0.99 


-0.034 


Illfe 


200 


-420 


-260 


-17000 


-11000 


4.0 


-0.034 



Effective parameters 
Set uo uj g gz 

la 0.0 1.0 0.5 -0.18 

Ha 0.0 1.0 1.0 -0.18 

nia 0.0 1.0 2.0 -0.18 

lb 0.0 1.0 0.5 -0.50 

life 0.0 1.0 1.0 -0.50 

nib 0.0 1.0 2.0 -0.50 

TABLE I: Parameter sets for the data of Figs. [3] and H All 
numbers are in units of (27r)-MHz. 



A. Photon number and revival probability 

Assuming an atom-cavity coupling strength gc/ (2Tr) = 
200 MHz, and choosing a selection of different laser fre- 
quencies and intensities to realize effective parameters in 
the ultrastrong and deep strong coupling regime, we plot 
in Fig. [3] the time evolution of the probability of being 
in the initial state, chosen to be |eO) = 1) <E) |0)cav, 
and the mean intracavity photon number ^a^a^. The 
physical parameters used and the effective parameters 
they give rise to are shown in Table HI To consider the 
deviation from the Rabi model Hamiltonian due to the 
(72-term in Eq. (|21l) . we compare two different values: 
(7z/(27r) = —0.18 and —0.5 MHz. We also compare two 
different field decay rates, k/{2tt) — 0.1 and 0.01 MHz, 
while for the full model, Eq. the spontaneous decay 
rate is 7/(27r) = 5.7 MHz. 

Over the timescale shown in Fig. [31 the effects of 
atomic spontaneous emission are negligible, and the ef- 
fective model fits the full model more or less perfectly 
with these parameters. We further see that, for the cho- 
sen values of k, the Hamiltonian dynamics dominate and 
there are clear signatures of the deep strong coupling 
regime, i.e., characteristic revivals of the initial state and 
(periodically) large intracavity photon numbers [20j . For 
comparison, we have also included in Fig.[3]the evolution 
of Eq. with gz = 0, i.e., the idealized (but damped) 
Rabi model system, shown as dotted lines. 

For longer evolution times, atomic spontaneous emis- 
sion must necessarily have some effect, both within the 
reduced system state space {|e), \g)} and by causing 
"leakage" into the other physical ground states (e.g., 
|2, — 1) or |1, 0)) as well. To consider this effect in the an- 
alytical model, one can include the spontaneous emission 
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5z/w = -0.18 5,/cj = -0.5 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 
t (/is) t {fis) 



FIG. 3: (Color online). Time evolution of the probability 
of the system being in the initial state, |eO), and the mean 
intracavity photon number (a^a). The two top rows are for 
the parameter sets denoted I in Tabled the two middle rows 
for the sets II, and the two bottom rows for the sets III. 
The left column corresponds to sets denoted by a, and the 
right column to b. The green solid lines are the solutions of 
Eq. ([9]), and the dashed black lines are those of the effective 
model Eq. (|22p . The two solutions coincide perfectly to the 
precision of the figure. The dotted lines are for the effective 
model with the same parameters, but with set to zero. 
Each panel shows results for two different h/{2tt) = 0.1,0.01 
MHz - the higher curves correspond to smaller k. 



term of Eq. ^ when doing the adiabatic elimination. 
This leads to additional dissipative terms (of Lindblad 
form) that connect the different physical ground states 
and are of order O ((7/2)if'^/A^) , where A is one of 
{Ai_2, Ai_2 ± W2}. We do not present these terms here, 
but rather investigate effects of spontaneous emission nu- 
merically in the full model. In particular, for parameter 



sets I-IIIa in Table HI we plot in Fig. |4] the long-time evo- 
lution of the total probability for the atom to be in \g) 
or |e), given an initial system state |eO). The decay rates 
of this probability are found to be at least several orders 
of magnitude smaller than the rates characterizing the 
effective Rabi model dynamics. 
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FIG. 4: (Color online). Population decay of the {|<7),|e)} 
atomic subspace in the simulation of the full system, with 
7/ {2tt) = 5.7 MHz and parameter sets la (green solid line), Ila 
(blue dashed line), and Ilia (black dotted line). Exponential 
fits of the form c-l-aexp(-n) yield F = 4.7 x 10~^ MHz (la), 
1.9 X 10"* MHz (Ha), and 7.8 x 10"* MHz (llla). Note that 
F scales linearly with (i — 1,2). 

While the cavity coupling constant in the above sim- 
ulations is well within projected values for future micro- 
toroidal resonators [2^, it is interesting to consider still 
smaller values in order to explore the requirements on gc 
for the simulation to be effective. We see from Eqs. ([M])- 
([TS]) that a smaller means the ratio of laser intensities 
to detunings must be larger in order to achieve the same 
effective coupling strength g. This means, in general, 
that the adiabatic elimination is less likely to be valid 
and spontaneous emission will be more significant. It 
follows that, for any gc, there is a limit to how far we can 
successfully push the effective g, or in other words the 
simulation scheme, into the ultrastrong regime. 

In Fig. Owe plot the time evolution of the mean pho- 
ton number and the revival probability of the initial state 
for a physical coupling constant gc/(27r) = 50 MHz and 
parameter sets as detailed in Table HIl When comparing 
with the previous parameter sets given in Table HJ we 
have essentially taken sets l-llla, but halved the detun- 
ings A 1.2 and roughly doubled the laser intensities rii,2, 
so that the effective coupling strengths remain the same: 
g/{2-K) = 0.5, I.O, and 2.0 MHz respectively. The param- 
eters 5 and 5c are chosen such that ~ 0.0 MHz and 
uj/{2tt) = 1.0 MHz, as before, while the value for gz comes 
out to gz/{'2.Ti) — —0.04 MHz (with this value of gz the 
solution of the effective model is indistinguishable from 
the case gz = 0). Fig.[5]shows that the simulation scheme 
still performs very well, at least on the timescale shown, 
for g/uj — 0.5 and g/uj = 1.0, which are both well into 
the ultrastrong coupling regime. For g/uj = 2.0, how- 
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Cavity and laser parameters 



Set Qc 0,1 Ai 



S 5c 



/lu = -0.04 



I 50 -390 -230 -16000 -10000 4.4 0.93 

II 50 -784 -470 -16000 -10000 18 0.93 

III 50 -1560 -940 -16000 -10000 70 0.93 

TABLE II: The parameter sets for the data of Figure [S] All 
numbers are in units of (27r)-MHz. These sets give effective 
parameters loo/{2tv) = MHz, tj/(27r) = 1 MHz, 5z/(27r) = 
-0.04 MHz, and 3/(27r) = 0.5 MHz (I), 1.0 MHz (II), 2.0 MHz 
(HI). 



ever, it is clear that the effective model breaks down. In 
particular, we find a substantial loss of population from 
the atomic subspace {|e), \g)} as a result of spontaneous 
emission. 



B. Field-atom entanglement 

We conclude this section by looking at the entangle- 
ment between the two level system {|g), |e)} and the field 
mode. Due to the realization of the Rabi model in the ul- 
trastr ong coupling regime in recent experiments in circuit 
QED |10l - [l3| , there has been interest in its application to 
potential quantum information technologies. In contrast 
to what one would expect from applying a rotating wave 
approximation (i.e., using Eq. ([2])) the ground state of 
the Rabi model is in general entangled |23| . and more so 
with stronger coupling g, as one might intuitively expect. 
It is therefore interesting to consider to what extent our 
effective model can be used to produce entangled states 
between atom and field. The dissipative model, Eq. , 
does not, however, evolve the system to the ground state 
of ifoff (see [23,[H,[32| for discussion and effective dissipa- 
tive models that do project onto the Rabi model ground 
state), and the mixing due to the dissipation limits the 
entanglement in the system, particularly in the long-time 
limit. 

As a quantitative measure of atom-field entanglement 
we use the logarithmic negativity, En{p) jicj. [ilf. which 
gives us an upper bound on the amount of distillablc 
entanglement present. In Fig. [5] we plot the time evo- 
lution of Epf{p), starting from the unentangled initial 
state |g0), for = 0.0, a;/(27r) = 1.0 MHz, .9j(27r) = 
-0.18 MHz, g/(27r) = {0.5,1.0,2.0} MHz and k/{2tt) = 
{0.1,0.01,0.0} MHz. After a short time proportional 
to g"^, the system evolves to a highly entangled state 
closely approximating the form jV') = {\e){\a) + \— a)) — 
\g)(\a) — I — q;))}/2, where \a) denotes a coherent state 
of the field mode of amplitude a (for g/{2TT) = 2.0 MHz 
with K/(27r) = 0.1 MHz, a fidelity of 0.98 is achieved). If 
the lasers were to be turned off at this time and the cav- 
ity field allowed to decay, then entanglement will persist, 
but between the atom and the light pulse propagating 
in the cavity output field. Alternatively, if the lasers 
were turned off and the atomic state measured (by, e.g.. 




t ins) 



FIG. 5: (Color online). Time evolution of the probability 
of the system being in the initial state, |eO), and the mean 
intracavity photon number (a^a). The two top panels are for 
the parameter set denoted I in Table [III the two middle panels 
for the set II, and the two bottom panels for the set HI. The 
green solid lines are the solutions of Eq. @, and the dashed 
black lines are those of the effective model Eq. (|22|) . Each 
panel shows results for two different K/{2n) = 0.1,0.01 MHz 
- the higher curves correspond to smaller k. 



fluorescence detection) , then the cavity mode will be pro- 
jected into one or other of the "Schrodinger cat" states 
\^P±) (x\a)±\-a) [281]. 

With K > the logarithmic negativity eventually de- 
cays towards zero for longer evolution times at a rate that 
actually grows with g. While a larger g gives a higher 
maximum entanglement at short times, it also makes 
the system evolution more sensitive to dissipation. This 
can be related to the state of the field mode pcav, which 
we examine in Fig. [7] with snapshots of the time evolu- 
tion of the Wigncr function for the cavity field for the 
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FIG. 6: Evolution of the logarithmic negativity, En{p), for 
the effective model Eq. (|22|l with initial state \gO). We con- 
sider three different coupling strengths: g/{2n) = 2.0 (solid 
lines), f.O (dashed lines) and 0.5 (dotted lines) MHz. Top 
panel: K/(27r) = O.f MHz, middle panel: k/{2tt) = O.Of MHz, 
bottom panel: k/{2tt) = 0.0 MHz. Other parameters are 
LOO = 0.0, uj/{2-k) = f.O MHz, and g,/{2Tv) = -O.fS MHz. 
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FIG. 7: (Color online). Contour plots of the Wigner function 
W{a), with a = (x + iy)/V2, showing snapshots of the cavity 
field's time evolution. g/(27r) = 2.0 MHz, K/i2n) = O.f Mhz, 
and other parameters as given in the text. Starting from the 
top left and traversing rows first, the times for the snapshots 
are t = {0, 0.25, 0.5, 0.75, f.O} /is. The corresponding mean 
photon numbers are (Ja^a) = {0, 6.8, f f .8, 5.5, 0.9}. 



{g, k}/{2tt) = {2.0, 0.1} MHz case. The Wigner function 
is defined by 

Wia) = -tT[D\a)p,,,D{a)i-lf-l (23) 

TT 

where D{a) = exp(aa^ — a*a) is the dispacenient opera- 
tor and pcav is the reduced density operator for the cav- 
ity field. Starting from the vacuum, the Wigner function 
evolves into two well separated peaks, attaining a max- 
imum separation at t = 0.5 /is, before the peaks nearly 
recombine close to the vacuum state after a full period, 
t — 1.0 ^s. The larger the coupling g, the greater sep- 
aration of the field amplitudes corresponding to the two 
peaks of the Wigner function and the greater sensitiv- 
ity of the coherence between these distinct field states to 
cavity decay. 

IV. THE NON-LINEAR ATOM-PHOTON 
INTERACTION 

In this section we continue our investigation of the ef- 
fective model, Eq. (|22p . but now examine in more detail 
the influence of the effective parameter on the sys- 
tem's behavior. The non-linear atom-photon interaction. 



g^cTza'' a/2, can be thought of as giving an effective, dy- 
namic shift to the cavity frequency, tu — >■ w = uj + gz(7z/2, 
or, alternatively, to the atomic frequency, — >■ wo = 
1^0 + gzO-^o,. For small gz, the dynamics is qualitatively 
similar to the Rabi model, as we have seen, but as \gz\ 
grows in magnitude towards the value 2^ an instability 
develops in the system that one can associate, at least 
qualitatively, with the potential for the effective cavity 
frequency a) to approach zero. In fact, if one considers 
just the generalised Hamiltonian, Hos, then for j g^j > 2a; 
the theoretical ground state is one of infinite photon num- 
ber and atomic state |e) {gz < —2a;) or \g) {gz > 2a;). 

Dissipation fundamentally alters this (unphysical) pic- 
ture, but signatures remain in the form of dramatic 
changes in the system properties around the "critical 
points" gz = ±2a; and distinct "phases" . In Fig. [5] we 
plot the steady state atomic inversion (cr^), intracavity 
photon number (^a^a), and intensity correlation func- 
tion .g'^''(0) = {a^ aa) / {a^ a)^ as a function of gz- The 
other parameters used are uj/{2'k) = a;o/(27r) = 1.0 MHz, 
k/{2tt) = 0.2 MHz, and g/{2TT) = {0.5, 1.0, 2.0} MHz. 

For sufficiently large the behavior around \gz\ = 2uj 
of each of the quantities plotted shows rapid and pro- 
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FIG. 8: Steady state atomic inversion, photon number and 
intensity correlation function g^^\Q) as a function of for 
p/(27r) = 0.5 MHz (dotted lines), 1.0 MHz (dashed lines) 
and 2.0 MHz (solid lines). Other parameters are u}/{2ti) = 
tJo/(27r) = 1.0 MHz and k,/{2'k) = 0.2 MHz 



nounccd changes. This critical- type behavior is further 
exemplified by the Wigner function of the cavity field, 
which we plot in Fig. ^ for g/(27r) = 2.0 MHz and a 
series of values of g^. In particular, over a very small 
range of gz values close to g^ = — 2w, W{a) changes 
from a weakly-split doublet aligned along the j/-axis to a 
strongly-resolved doublet aligned more closely with the x- 
axis, with an intermediate ("co-existent") phase in which 
VF(a) exhibits four distinct peaks. With a finite value of 
Wo, the properties of the system are not symmetric in gz 
and, in particular, the Wigner function displays only a 
single maximum in the region gz > 2lu. 

Interestingly, a semiclassical model of the system also 
reveals distinct phases about critical values gz — ±2w. 
In fact, for the region \gz\ < 2uj the semiclassical treat- 
ment gives results in good quantitative agreement with 
the results presented above, while for gz < —2a; at least 
qualitative agreement is found. For gz > 2lu it would ap- 
pear that no stable semiclassical solutions exist, but we 
defer a detailed investigation of the semiclassical dynam- 
ics, and its comparison with quantum results for both the 
single- and multiple-atom system, to a future work. 

The value of gz also strongly influences the time- 
dependent behavior of the system. In Fig.[TU]we plot the 
atomic inversion and cavity photon number as a function 



Si 



Si 




FIG. 9: (Color online). Contour plots of the Wigner func- 
tion W(a), with a — (x + iy)l\f2, for varying values of g^, 
with iJo/(27r) = tj/(27r) = 1.0 MHz, ff/(27r) = 2.0 MHz, and 
K/(27r) = 0.2 MHz. 



of time, when starting from an initial state |e, 0), for five 
different values of gz. We compare the effective model, 
Eq. dill), with the full *^Rb model, Eq. dH), and find ex- 
cellent agreement over the timescale shown. The physical 
and effective parameters used are given in Table lllll The 
figure shows, in particular, that characteristic timescales 
for oscillations and decay towards a steady state vary 
significantly with as one expects for a system exhibit- 
ing critical-type behavior and manifestly distinct phases 
associated with variations of this parameter. 
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Cavity and laser parameters 



Set Qc ^1 ^2 Ai 



5^ 



Iz/lo = —3.0 



g^/uj = -3.0 



I 200 -120 -40 -9700 -3700 1.4 -1.9 -3.0 

II 200 -130 -53 -11000 -4800 1.2 -1.2 -1.5 

III 200 -320 -240 -26000 -20000 1.8 0.40 -0.18 

IV 200 130 53 -11000 -4800 -0.80 -3.2 1.5 

V 200 120 40 -9700 -3700 -0.65 -3.9 3.0 

TABLE III: The parameter sets referred to in Section |IV] All 
numbers are in units of (27r)-MHz. The effective parameters 
realized are in each case uo/{2n) = u}/{2-k) = g/{2n) = 1.0 
MHz, while is varied as given in the rightmost column. 
Notice that for set IV and V, we have used Hcs — —//off in 
Eq. 



V. CONCLUSION 

We have proposed a scheme for realizing a system 
with dynamics described by an effective generalized Rabi 
model. This scheme is based on resonant Raman transi- 
tions between stable ground states of a ^^Rb atom. By 
numerical solution of the master equation, including cav- 
ity decay and spontaneous emission, we have verified the 
feasibility of implementing this scheme, given sufficiently 
strong coupling of the atom to a high finesse optical 
cavity mode, as should be achievable with microtoroidal 
whispering gallery mode resonators. This scheme offers 
a means of performing a quantum simulation of the Rabi 
model with coupling strengths in both the ultrastrong 
[g < uj) and deep strong {g > w) regimes, which is inter- 
esting both from a fundamental point of view and in the 
context of preparing highly nonclassical states. 

Our scheme also offers a generalization of the Rabi 
model, with the addition of a non-linear atom-photon 
coupling in the Hamiltonian that can be interpreted as 
a dynamical shift of the cavity mode frequency. This 
coupling, quantified by the parameter (7^, leads to funda- 
mentally new phases of system behavior for magnitudes 
of gz exceeding twice the effective cavity mode frequency. 
Transitions between these phases are especially sharp in 
the deep strong coupling regime, with clear signatures 
in the properties of the cavity output field (i.e., in the 
output intensity and intensity correlation function). Ex- 
tension of these results to the multiple-atom case should 
offer new opportunities for the study of critical phenom- 
ena in many-body cavity QED. 

Finally, a possibility that we have not explored in this 
paper, but would like to point out, is that the effective 
parameters in our model could easily be made time de- 
pendent. A modulation of the effective coupling constant 
p, for example, could be introduced by varying the laser 
Rabi frequencies r2i_2 (i.e., by varying the laser intensi- 
ties). In the context of the Rabi model, this has been 
predicted to have some fascinating consequences, such as 
the generation of photons from the vacuum p^ . In fact, 
this can be viewed as being an analog of the dynamical 
Casimir effect, and it would be interesting to consider 




FIG. 10; (Color online). Time evolution of the atomic inver- 
sion, {(T2) (solid green lines in left column), and the intracavity 
photon number (a) a}f (solid green lines in right column) for 
Eq. Q. The dashed lines are the corresponding results for 
the effective model, Eq. (|22|) . The solutions are nearly indis- 
tinguishable. The effective parameters are in each case set to 
ujo/(2-k) = uj/{2-k) = g/(27r) = 1.0 MHz, k,/(2-k) = 0.2 MHz, 
7/(27r) = 5.7 MHz, while is varied. The corresponding 
cavity and laser parameters used are given in Table Hill with 
set I to V starting from the top to bottom row. 



how our scheme could be used to simulate this elusive 
phenomenon. 
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